Generalizing the autonomous Kepler Ermakov system in a 

Riemannian space 



(N 

o 

(N 



c3 



X 



Michael Tsamparlis* and Andronikos Paliathanasis^ 

Faculty of Physics, Department of Astrophysics - Astronomy - Mechanics, 
University of Athens, Panepistemiopolis, Athens 157 83, GREECE 

Abstract 



We generalize the two dimensional autonomous Hamiltonian Kepler Ermakov dynamical system to three 
dimensions using the sl(2,R) invariance of Noether symmetries and determine all three dimensional au- 
l/") | tonomous Hamiltonian Kepler Ermakov dynamical systems which are Liouville integrable via Noether sym- 

metries. Subsequently we generalize the autonomous Kepler Ermakov system in a Riemannian space which 
admits a gradient homothetic vector by the requirements (a) that it admits a first integral (the Riemannian 
■ Ermakov invariant) and (b) it has sl(2, R) invariance. We consider both the non-Hamiltonian and the Hamil- 

Oh| tonian systems. In each case we compute the Riemannian Ermakov invariant and the equations defining the 

(— ( ■ dynamical system. We apply the results in General Relativity and determine the autonomous Hamiltonian 

Riemannian Kepler Ermakov system in the spatially flat Friedman Robertson Walker spacetime. We con- 
sider a locally rotational symmetric (LRS) spacetime of class A and discuss two cosmological models. The 
first cosmological model consists of a scalar field with exponential potential and a perfect fluid with a stiff 
equation of state. The second cosmological model is the f(R) modified gravity model of A{, C CDM. It is shown 
CN) , that in both applications the gravitational field equations reduce to those of the generalized autonomous 

Riemannian Kepler Ermakov dynamical system which is Liouville integrable via Noether integrals. 
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1 Introduction 



The Ermakov system has its roots in the study of the one dimensional time dependent harmonic oscillator 

X + UJ 2 (t)x = 0. (1) 

Ermakov .1] obtained a first integral J of this equation by introducing the auxiliary equation 

p + uj 2 (t)p = p~ 3 (2) 
eliminating the 0J 2 (t) term and multiplying with the integrating factor px — px : 

J=\[{px-pxf + {x/pf]. (3) 

The Ermakov system was rediscovered nearly a century after its introduction [2] and subsequently was 
generalized beyond the harmonic oscillator to a two dimensional dynamical system which admits a first integral 
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[3]. In a series of papers the Lie, the Noether and the dynamical symmetries of this generalized system have 
been studied. A short review of these studies and a detailed list of relevant references can be found in [4]. 
Earlier reviews of the Ermakov system and its numerous applications in divertive areas of Physics can be found 
in 

The general Ermakov system does not admit Lie point symmetries. The form of the most general Ermakov 
system which admits Lie point symmetries has been determined in [7] and it is called the Kepler Ermakov 
system [3 [9]. It is well known that these Lie symmetries are a representation of the sl(2,R) algebra. 

In an attempt to generalize the Kepler Ermakov system to higher dimensions Leach [5] used a transformation 
to remove the time dependent frequency term and then demanded that the autonomous 'generalized' Kepler 
Ermakov system will posses two properties (a) A first integral, the Ermakov invariant and (b) sl(2, R) invariance 
to Lie symmetries. It has been shown that the invariance group of the Ermakov invariant is reacher than 
sl(2,R) [TU]. The purpose of the present work is to use Leach's proposal and generalize the autonomous 
Kepler Ermakov system in two directions (a) to higher dimensions using the sl(2,R) invariance with respect 
to Noether symmetries (provided the system is Hamiltonian) and (b) in a Riemannian space which admits a 
gradient homothetic vector (HV). 

The generalization of the autonomous Kepler Ermakov system to three dimensions using Lie symmetries has 
been done in [5] . In the following we use the results of [TT] to generalize the subset of autonomous Hamiltonian 
Kepler Ermakov systems to three dimensions via Noether symmetries. We show that there is a family of three 
dimensional autonomous Hamiltonian Kepler Ermakov systems parametrized by an arbitrary function / which 
admits the elements of sl{2, R) as Noether symmetries. Each member of this family admits two first integrals, 
the Hamiltonian and the Ermakov invariant. 

We use this result in order to determine all three dimensional Hamiltonian Kepler Ermakov systems which 
are Liouville integrable via Noether symmetries. To do this we need to determine all members of the family, 
that is those functions /, for which the corresponding system admits an additional Noether symmetry. 

The results of III indicate that there are two cases to be considered i.e. Noether point symmetries resulting 
from linear combinations of (a) translations and (b) rotations (elements of the so (3) algebra). In each case we 
determine the functions / and the required extra time independent first integral. 

The above scenario can be generalized to an n- dimensional Euclidian space as Leach indicates in [5] , however 
at the cost of major complexity and number of cases to be considered. Indeed as it can be seen by the results 
of the situation is complex enough even for the three dimensional case . 

We continue with the generalization of the Kepler Ermakov system in a different and more drastic direction. 
We note that the Ermakov systems considered so far are based on the Euclidian space, therefore we may call them 
Euclidian Ermakov systems. Furthermore the sl(2,R) symmetry algebra of the autonomous Kepler Ermakov 
system is generated by the trivial symmetry <9 t and the gradient HV of the Euclidian two dimensional space E 2 . 
Using this observation we generalize the autonomous Kepler Ermakov system (not necessarily Hamiltonian) in 
a n— dimensional Riemannian space which admits a gradient HV using either Lie or Noether symmetries. The 
new dynamical system we call the Riemannian Kepler Ermakov system. This generalization makes possible the 
application of the autonomous Kepler Ermakov system in General Relativity and in particular in Cosmology. 

Concerning General Relativity we determine the four dimensional autonomous Riemannian Kepler Ermakov 
system and the associated Riemannian Ermakov invariant in the spatially flat Freedman - Robertson - Walker 
(FRW) spacetime; we use the results of [EJUS] to calculate the extra Noether point symmetries. The applications 
to cosmology concern two models for dark energy on a locally rotational symmetric (LRS) space time. The 
first model involves a scalar field with an exponential potential minimally interacting with a perfect fluid with 
a stiff equation of state. The second cosmological model is the f(R) modified gravity model of h.b c CDM . It is 
shown that in both models the gravitational field equations define an autonomous Riemannian Kepler Ermakov 
system which is integrable via Noether integrals. 

The structure of the paper is as follows. In section [5] we review the main features of the two dimensional 
autonomous Euclidian Kepler Ermakov system. In section [3] we discuss the general scheme of generalization of 
the two dimensional autonomous Euclidian Kepler Ermakov system to higher dimensions and to a Riemannian 
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space which admits a gradient HV. In section [4] we consider the generalization to the 3d autonomous Euclidian 
Hamiltonian Kepler Ermakov system by Noether symmetries and determine all such systems which are Liouville 
integrable. In section [5] we define the autonomous Riemannian Kepler Ermakov system by the requirements 
that it will admit (a) a first integral (the Ermakov invariant) and (b) posses sl(2,R) invariance. In section 15.11 
we consider the non-conservative autonomous Riemannian Kepler Ermakov system and derive the Riemannian 
Ermakov invariant and in section 15.21 we repeat the same for the autonomous Hamiltonian Riemannian Kepler 
Ermakov system. In the remaining sections we discuss the applications of the autonomous Hamiltonian Rie- 
mannian Kepler Ermakov system in General Relativity and in Cosmology. Finally in section [5] we draw our 
conclusions. 



2 The two-dimensional autonomous Kepler Ermakov system 

In [7] Hass and Goedert considered the most general 2d Newtonian Ermakov system to be defined by the 
equations: 

x + uj 2 (t 7 x 7 y,x,y)x = -^f (-) (4) 
yx z \xJ 

y + u 2 (t,x,y,x,y)y = — (~) ■ (5) 
xy z \xJ 

This system admits the Ermakov first integral 

i rv/ x rv/x 

I=-( xy -yx) 2 + f(r)dr+ g(r)dr. (6) 



If one considers the transformation: 



n 2 = lu 2 - 

xy 



^(") 

y 3 \x/ 



H-)=/(-)-^(- 



then equations (HJ),© take the form: 



x + il 2 (x,y,x,y)x=^-F(^-) (7) 
x z y \xJ 

y + Q 2 (x,y,x,y)y = 0. (8) 

Due to the second equation, except for special cases, the new function 51 is independent of t and depends 
only on the dynamical variables x, y and possibly on their derivative. The Ermakov first integral in the new 
variables is: ^ 

I=\{xy~yx) 2 + J"* F{\)d\. (9) 

The system of equations |7} , © defines the most general 2d Ermakov system and produces all its known forms 
for special choices of the function f2. For example the weak Kepler Ermakov system [5] defined by the equations 

i 

x + ^{t)x+^H{x,y)-\f( V -) = (10) 
r J x a \xs 

y + u 2 {t)y+\H{x lV )-^g( V -) = (11) 
r'' y a \x/ 



where H , /, g are arbitrary functions of their argument, is obtained by the function 

n 2 (x,y) =u 2 (t)+H(x,y)/r 3 (12) 
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with Ermakov first integral 

I = ~ (xy - yxf + J 1 [A/ (A) - X^g (A)] dX. (13) 

The weak Kepler Ermakov system does not admit Lie point symmetries, however the property of having a first 
integral prevails. The system of equations (ITU)) , (ITT1) admits the si (2,i?) as Lie point symmetries [13] only for 
h (-) 

H (x, y) — —p 2 r 3 H — *±A where p is either real or pure imaginary number. This is the Kepler Ermakov system 
defined by the equations 



»+(«*« - ° <»> 

It is well known (see [13]) that the oscillator term cj 2 (t) — /i 2 in (|14l) . ([T5j) is removed if one considers new 
variables T, X, Y defined by the relations: 

T = J p- 2 dt,X = p~ l x ,Y = p' x y (16) 

where p is any smooth solution of the time dependent oscillator equation 

p+(Lo 2 (t)-n 2 )p = 0. (17) 

In [13] it is commented that "the effect of ji 2 (C in the notation of [13]) is to shift the time - dependent 
frequency function". However this is true as long as u>(t) ^ 0. When u(t) = one has the autonomous Kepler 
Ermakov system whose Lie symmetries span the si(2, R) algebra with different representations for /i = and 

Before we justify the need for the consideration of the two cases p, = and p ^ 0, we note that by applying 
the transformation 

s= I v~ 2 dT , x = v~ 1 X , y = v~ 1 Y (18) 



^ • 4 = (19) 



where v satisfies the Ermakov Pinney equation 

dT 2 ' v 3 



to the transformed equations 



£ + iH£) - <20) 

d 2 Y 1 Y, fY\ 1 (Y 



dT 2 b?x \x J y 3 - y \X 

we retain the term p? and obtain the autonomous Kepler Ermakov system of 



5«h7 = o Pi) 



y - m 2 ^ 



+ 3^(9 -^(") - °- (23) 
r J a; \x/ y a \x/ 



The above transformations show that the consideration of the autonomous Kepler Ermakov system is not a 
real restriction. 

We discuss now the need for the consideration of the cases p = and p ^ 0. In [T^] we have determined the 
Lie symmetries of the autonomous 2d Kepler Ermakov system and we have found two cases. The first Case I 
concerns the autonomous Kepler Ermakov system with p — and has the Lie symmetry vectors 
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X= (c x + c 2 2t + c 3 t 2 ) d t + (c 2 + c 3 t) rd r (A» = 0) (24) 
The second Case II concerns the same system with /i ^ and has the Lie symmetry vectors 

Cl + c 2 -e 2 ^ - c 3 -e- 2 ^) d t + (c 2 e 2 ^ + c 3 e- 2 ^) rd r ( fi ? 0) (25) 

where in both cases rda = xd x + yd y is the gradient HV of the 2d Euclidian metric. Each set of vectors in 
(|25|) . ((24)) is a representation of the si (2, R) algebra and furthermore each set of vectors is constructed from the 
vector dt and the gradient HV rd r of the Euclidian two dimensional space E 2 . 

The essence of the difference between the two representations is best seen in the corresponding first integrals. 
If a Kepler Ermakov system is Hamiltonian then the Lie symmetries are also Noether point symmetries therefore 
in order to find these integrals we determine the Noether invariants. The Noether symmetries of the Kepler 
Ermakov system have been determined in |12] . For the convenience of the reader we repeat the relevant material. 

Equations (|22[). ([2"3")l follow from the Lagrangian [T3] 

where C(8) = c + sec 2 6 f (tan 8) + esc 2 8g (tan 9) provided the functions /, g satisfy the constraint: 

sin 2 Of (tan 6) + cos 2 8 g' (tan 8) = 0. (27) 
The Ermakov invariant in this case is |13j . 

J = r 4 8 2 + 2C (8) . (28) 
Because the system is autonomous the first Noether integral is the Hamiltonian [T^] 

E= 1 -(r 2 +r 2 8 2 )+ 1^ + ^(8) (29) 

In addition to the Hamiltonian there exist two additional time dependent Noether integrals as follows: 
\i = 

h = 2tE~rf (30) 
h = t 2 E - trr + K 2 (31) 

I[ = Q-E - rf + ^ e 2 ^ (32) 

I' 2 = Q-E + rf + fir 2 ^j eT 2 ^ . (33) 

We note that the Noether integrals corresponding to the representation (|24l) are linear in t whereas the ones 
corresponding to the representation (|25p are exponential. Therefore the consideration of the cases ji = and 
/i 7^ is not spurious otherwise we loose important information. This latter fact is best seen in the applications 
of Noether symmetries to field theories where the main core of the theory is the Lagrangian. In these cases the 
potential is given and, as it has been shown in \V2\ . a given potential admits certain Noether symmetries only; 
therefore one has to consider all possible cases. We shall come to this situation in section [7] where it will be 
found that the potential selects the representation (|25l) . 

To complete this section we mention that for a Hamiltonian Kepler Ermakov system the Ermakov invariant 
P5| is constructed [T3] from the Hamiltonian and the Noether invariants (|3"2"j) . (|3"3"|) as follows: 

J = E 2 - I[I' 2 . 

Finally in |14| it is shown that the Ermakov invariant is generated by a dynamical Noether symmetry of the 
Lagrangian (|26p . a result which is also confirmed in |15j . 
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3 Generalizing the autonomous Kepler Ermakov system 



We consider the generalization of the two dimensional autonomous Kepler Ermakov system [9j [3 [16j [TT1 [T8l [19] 
using a geometric point of view. From the results presented so far we have the following: 

(i) Equations Q , ([3]) which define the Ermakov system employ coordinates in the Euclidian two dimensional 
space, therefore the system is the Euclidian Ermakov system. 

(ii) The autonomous 2d Euclidian Kepler Ermakov system is defined by equations (|2"2"|) and (f2"3"|) 

(iii) The Lie symmetries of the Kepler Ermakov system span the sl(2, R) algebra. These symmetries are 
constructed from the vector <9 t and the gradient HV of the space E 2 . 

(iv) For the autonomous Hamiltonian Kepler Ermakov system the Lie symmetries reduce to Noether sym- 
metries and the Ermakov invariant follows from a combination of the resulting three Noether integrals, two of 
which are time dependent. Furthermore, the Ermakov invariant is the Noether integral of a dynamical Noether 
symmetry. 

The above observations imply that we may generalize the Kepler Ermakov system in two directions: 

a. Increase the number of dimensions by defining the n— dimensional Euclidian Kepler Ermakov system 
and/or 

b. Generalize the background Euclidian space to be a Riemannian space and obtain the Riemannian Kepler 
Ermakov system. 

Concerning the defining characteristics of the Kepler Ermakov system we distinguish three different prop- 
erties of reduced generality: The property of having a first integral; the property of admitting Lie/Noether 
symmetries and sl{2,R) invariance and the property of being Hamiltonian and admitting sl{2,R) invariance 
via Noether symmetries. 

Following Leach [9] we generalize the autonomous Kepler Ermakov system to higher dimensions by the 
requirement: 

The generalized autonomous (Euclidian) Kepler Ermakov system admits the sl(2,R) algebra as a Lie 

symmetry algebra. 

In order to exploit the significance of the sl(2, R) invariance we refer to two theorems, which relate the Lie 
and the Noether point symmetries of an autonomous dynamical system moving in a Riemannian space with the 
symmetries of the space (for a detailed statement of these theorems see |12j). 

Theorem 1 The Lie point symmetries of the equations of motion of an autonomous dynamical system moving 
under the force F 3 (x l ) in a Riemannian space with metric gij, namely 

x i + T i jk x j x k = F i (34) 

are given in terms of the generators of the special projective algebra of the metric gtj. 

If the equations of motion follow form the standard Lagrangian 

L {x 3 ,x J ) = ^g ij x i x j - V (x 3 ) (35) 

where V(x 3 ) is the potential, then the following theorem relates the Noether point symmetries of the Lagrangian 
with the symmetries of the metric. 

Theorem 2 The Noether point symmetries of the Lagrangian i35\) of an autonomous Hamiltonian system 
moving in a Riemannian space with metric g^ are generated from the homothetic algebra of g^ . 
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4 The 3d autonomous Euclidian Kepler Ermakov system 



The generalization of the autonomous Euclidian Kepler Ermakov system using si (2, R) invariance of Lie sym- 
metries has been done in [5J [TH1 HH] ■ In this section, using of the results of [TT] , we give the generalization 
of the autonomous Euclidian Hamiltonian Kepler Ermakov system to three dimensions by demanding sl(2,R) 
invariance with respect to Noether symmetries. The reason for attempting this generalization is that it leads 
to the potentials for which the corresponding extended systems are Liouville integrable. Furthermore indicates 
the path to the n— dimensional Riemannian Kepler Ermakov system. 

Depending on the value fi ^ or fi = we consider the three dimensional Hamiltonian Kepler Ermakov 
systems of type I and type II. 

4.1 The 3d autonomous Hamiltonian Kepler Ermakov system of type I (fx ^ 0) 

For /i ^ the admitted Noether symmetries are required to be (see ([25]) ): 

X 1 = dt, X± = -e ±2 ^d t ± e ±2tit Rd R . (36) 
M 

From table 6, line 2 for T = le* 2 ^' f [XT] we find that for these vectors the potential is V(R,4>,6) = 
— §--R 2 + -^f {6,6) hence the Lagrangian 

L=\(r 2 + R 2 tf + R 2 sin 2 6 9 2 ) + J^-R 2 ±f (0, 6) . (37) 
The equations of motion, that is, the equations defining the generalized dynamical system are 

R-R<j) 2 -R sin 2 6 6 2 -fj?R- 4r/ = ( 38 ) 

R d 

•• 2 • • -1 

R<p-sin6cos6 6 2 + — f t4> = (39) 



R r R 4 

2 • • • • 1 

— Rd + cot 6 66 H . 

R i? 4 sin 



Rd + cot 6 9<j> + — , 2 f >6 = 0. (40) 



The Noether integrals corresponding to the Noether vectors are 

1 ( hi i D2i2 i D 2 „j„2 x 02\ A* r>2 1 



E = -^ + R'6 2 +R'sm 2 6 9 2 )~^R z + —f(6,6) (41) 
1+ = -e 2 ^E - e 2 ^RR + fie^R 2 (42) 

J_ = -e-^E + e-^RR + fie-^R 2 (43) 
H 

where E is the Hamiltonian. The Noether integrals I± are time dependent. Following [T5] we define the time 
independent combined first integral 

J = E 2 - I+I- = R 4 <j) 2 + i? 4 sin 2 6 6 2 + 2f (9, 6) . (44) 

Using (|4"4"|) the equation of motion (|3"5|) becomes 

R-^R=± (45) 

which is the autonomous Ermakov- Pinney equation [20]. Therefore J is the Ermakov invariant [5]. 

An alternative way to construct the Ermakov invariant (1441) is to use dynamical Noether symmetries 22 . 
Indeed one can show that the Lagrangian (|3"T|) admits the dynamical Noether symmetry Xd — K^x^di where 
Kij is a Killing tensor of the second rank whose non- vanishing components are — i? 4 , Kgg = R 4 sin 6. The 
dynamical Noether symmetry vector is Xd = R 2 {^d^ + Odg^j with gauge function 2/ (6, 6). 
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4.2 The 3d autonomous Hamiltonian Kepler Ermakov system of type II (// = 0) 

For [i = the Noether symmetries are required to be (see (|2"4"|) ) [3]: 

X 1 = d t , X 2 = 2td t + Rd R , X 3 = t 2 d t + tRd R . (46) 

From table 5, line 3 for T(t) = 1 and table 6, line 1 for T(t) =toi 11 we find that the potential V (R, <j>, 9) = 
■jjsf (9, (ft) hence the Lagrangian 

L' = \ (it 2 + R 2 <P 2 + R 2 sin 2 ft9 2 ) - ^ / (6, 4>) . (47) 
The equations of motion are (j3"8")l - (|4U1) with /i — 0. The Noether invariants of the Lagrangian (|4"Tf are 

E = ^[r 2 +R 2 4> 2 +R 2 s\n 2 <i ) e 2 ) +-Lf(0,4>) (48) 
h = 2tE' - RR (49) 
h = t 2 E' -tRR+^R 2 . (50) 

We note that the time dependent first integrals ii ; 2 are linear in t whereas the corresponding integrals I± 
of the case fj, ^ are exponential. From Ii^ we define the time independent first integral J = AI2E' — I 2 which 
is calculated to be 

J = R A cj) 2 + R 4 sin 2 cf> 8 2 + 2f (9, <p) . (51) 

Using (|5Tj) the equation of motion for R (t) becomes R — -jm = which is the one dimensional Ermakov- 
Pinncy equation, hence J' is the Ermakov invariant [9]- As it was the case with the three dimensional 
Hamiltonian Kepler Ermakov system of type I, the Lagrangian (|47[) admits the dynamical Noether symme- 

try X D = R 2 (0S + 9d e ) whose integral is the ((ST]) . 

4.3 Integrability of the 3d autonomous Hamiltonian Euclidian Kepler Ermakov 
system 

The 3d autonomous Hamiltonian Euclidian Kepler Ermakov systems need three independent first integrals in 
involution in order to be Liouville integrable. As we have shown each system has the two Noether integrals 
(E, J) , hence we need one more Noether symmetry. Such a symmetry exists only for special forms of the 
arbitrary function / (9, (ft). From tables 5, 6 and 9-11 of [UJ we find that extra Noether symmetries are possible 

onlj0 for linear combinations of translations (i.e. vectors of the form aA c>A where a A are constants) and/or 

A=l 

rotations (i.e. elements of so(3)). 

4.3.1 Noether symmetries generated from the translation group 

We determine the functions / (9, (ft) for which the 3d autonomous Hamiltonian Euclidian Kepler Ermakov 
systems admit extra Noether symmetries for linear combinations of the translation group. 
The Lagrangian (|3"7) 

In Cartesian coordinates the Lagrangian (l37l) is 

L (.H..H) = 1 -(i 2 +y 2 + z 2 ) + £ (* 2 + v 2 + z 2 ) - 1/, (|, J) (52) 

1 The linear combination of an element of so(3) with a translation does not give a potential, hence an additional Noether 
symmetry. 



where fi — (l + ^ + . From table 6, line 1 of [11 with m = —fi 2 , p = we find that the Lagrangian 

(|52p admits Noether symmetries which are produced from a linear combination of translations, if the function 
fx (-> -) has the form 

*(S-5)-^V(^)- 

In this case the Lagrangian (15"2"j) admits at least the following two extra Noether symmetries (see table 9, 
line 1 of [IT]) 

3 

*± =e if "^a A 9A (54) 

A=l 

with corresponding Noether integrals 

I± = (iX^J T ^ (E oA ^J • ( 55 ) 

We note that the first integrals I± are time dependent; however the first integral 

J 2 = I+I- = (ax + by + cz) 2 + y 2 (ax + by + cz) 2 (56) 



is time independent. As it was the case with the Ermakov invariant (1441) the integral J 2 is possible to be 
constructed directly from the dynamical Noether symmetry X' D = Kfy ji l di, where K( 2 )ij is a Killing tensor 
of the second rank [22, 23 , with non- vanishing components 

K n = a 2 , K 22 = b 2 , K 33 = c 2 
K {12] = 2ab , K (13} = 2ac , K {23) = 2bc 

so that the dynamical symmetry vector is 

X' D = (a 2 + ab + acj xd x + (b 2 + ab + be) yd y + (c 2 + ac + 6c) zd z . (57) 

The Ermakov invariant J (see (|44p ) in Cartesian coordinates is 

J = 2E (x 2 + y 2 + z 2 ) - (xx + yy + zz) 2 . (58) 

The first integrals J, J 2 are not in involution. Using the Poisson brackets we construct new first integrals 
and at some stage one of them will be in involution. These new first integrals can also be constructed form 
corresponding dynamical Noether symmetries. 

An example of a known Lagrangian of the form (fB"2")) is the three body Calogero-Moser Lagrangian [2"Tl |2~4"1 12"5] 

L = \(x 2 + y 2 + z 2 )-^(x 2 + y 2 : + ^- T ± -2--^-- ± ^- (59) 
1 1 (x-y) (x-z) (y-z) 

The extra Noether symmetries of this Lagrangian are produced by the vector (|5"4")l for a A — (1,1,1). 
The Lagrangian (|47| 



In Cartesian coordinates the Lagrangian (|47p is 

L (xt,*>) - \ (x 2 + y 2 + z 2 ) \f H ( V ~, Z -\ . (60) 
v ' 2 X ' x l \x xJ 

According to table 10, line 1 and table 11, line 1 of [TT] (with m = , p = 0) the Lagrangian (jBT)]) admits extra 
Noether point symmetries for a linear combination of translations if the function / is of the form (l53l) . In this 
case the corresponding Noether integrals are 

3 33 
/( = ^a A x A , I' 2 = t^2a A x A - ^a A x A - (61) 

A=l A=l A=l 
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Example of such a Lagrangian is the Calogero-Moser Lagrangian [21] (without the oscillator term) 
L = \(x 2 + y 2 + z 2 ) — j — 2 " — 2 • 

2 (x-y) (x-z) (y-z) 



(62) 



For the Lagrangian (jf)2"]) we have the first integrals E, J, /j- From the integrals J, I[ we construct the first 
integral $ = , {J, /{}} . It is easy to show that the integrals E,I[,$ are in involution hence the dynamical 
system is Liouville integrable. We remark that the first integrals E, J, ij,-^ can & l so be computed by making 
use of the Lax pair tensor [23] . 



4.3.2 Noether symmetries generated from so (3) 

The elements of so (3) in spherical coordinates are the three vectors CK 1 ' 2 ' 3 

CK 1 = sin 8drh + cos 9 cot d>dg , CK 2 = cos Odd, — sin 9 cot d 



CK = d ff 



(63) 



which are also the KVs of the Euclidian sphere. 

In this case the symmetry condition c.f. (eq.(56) of 



^CK 



1 



f ( 



or, equivalently 



1 



p = Q 

where gty is the metric of the Euclidian sphere, that is 

ds 2 = d<h 2 + 



) becomes 
-p = 

n 2 d> d9 2 . 



(64) 



(65) 



(66) 



We infer that the problem of determining the extra Noether point symmetries of Lagrangian (|37p generated 
from elements of the so (3) is equivalent to the determination of the Noether point symmetries for motion on 
the 2d sphere, a problem which has been answered in [TZ] . 

It is easy to show that the integrals of table 7 of [11] are in involution with the Hamiltonian and the Ermakov 
invariant, hence the system is Liouville integrable via Noether symmetries. 

The above results are extended to the case in which the system moves on the hyperbolic sphere that is, it 
has Lagrangian 



R 



2 <P 2 + R 2 sinh 2 6 9 2 )+ ^R 2 - -L 



(67) 



We reach at the following conclusion. 



Proposition 3 The three dimensional autonomous Hamiltonian Kepler Ermakov system with Lagrangian \37\ l 
is Liouville Integrable via Noether point symmetries which are generated from a linear combination of the three 
elements of the so (3) algebra, if and only if the equivalent dynamical system in the fundamental hyperquadric^ 
of the three dimensional flat space is integrable. 

We note that it is possible a three dimensional autonomous Kepler Ermakov system to admit more Noether 
symmetries, which are due to the rotation group and the translation group (but not to a linear combination of 
elements from the two groups). For example the 3d Kepler Ermakov system with Lagrangian [27] 



L = Ui 2 +y 2 



z 2 )- 



2 See 1261 for definitions. 
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has the following extra Noether point symmetries (in addition to the elements of si (2,R)) 

y 1 = d x + d v , Y 2 = d x + d z 



Y 3 = t(d X +dy), Y 4 =t(d X + d Z ) 



Y" = (y- z )d x -(x + z)d y + (x + y)d z . 

The vectors Y 1,2 and Y 3 ' 4 follow from ((54|) for a A — (1, 1, 0) and a A — (1, 0, 1) respectively whereas Y 5 is a 
linear combination of the three elements of so (3). The Noether integrals of the Noether symmetries V 1-5 are 
respectively 

lY\ = x + y , Iy 2 = x + z (69) 
I Ys = t (x + y) - (x + y) , I Yi = t (x + z ) - (x + z) (70) 
Iy 5 = (y-z)x-(x + z)y + (x + y)z. (71) 

It is clear that in order to extend the Kepler Ermakov system to higher dimensions one needs to have the type 
of results of [TT]; therefore the remark made in [S] that the 'notion is easily generalized to higher dimensions' 
has to be understood as referring to the general scenario and not to the actual work. 

5 The autonomous Riemannian Kepler Ermakov system 

As it has been noted in section [5] the Kepler Ermakov systems considered so far in the literature are the New- 
tonian Kepler Ermakov systems. In this section we make a drastic step forward and introduce the autonomous 
Riemannian Kepler Ermakov system of dimension n. The generalization we consider is based on the following 
definition 

Definition 4 The n— dimensional autonomous Riemannian Kepler Ermakov system is an autonomous dynam- 
ical system which: 

a. R is defined on a Riemannian space which admits a gradient HV 

b. Admits a first integral, which we name the Riemannian Ermakov first integral and it is characterized by the 
requirement that the corresponding equation of motion takes the form of the Ermakov Pinney equation. 

c. R is invariant at least under the sl{2, R) algebra, which is generated by the vector dt and the gradient HV of 
the space. 

There are two types of n— dimensional autonomous Riemannian Kepler Ermakov systems. The ones which 
are not Hamiltonian and admit the sl(2, R) algebra as Lie symmetries and the ones which are Hamiltonian and 
admit the sl(2,R) algebra as Noether symmetries. 

5.1 The non Hamiltonian autonomous Riemannian Kepler Ermakov system 

Consider a n— dimensional Riemannian space which admits a gradient HV. It is well known that the metric of 
this space can always be written in the form [23 [29] 



ds 2 = du 2 + u 2 h AB dy A dy B (72) 
where the Latin capital indices A, B, .. take the values 1, . . . , n — 1 and Hab — ^ab (y C ) is the generic n — 1 

b 



metric. The gradient HV of the metric is the vector H l = ud u (ip — 1) generated from the function H = \u 2 



We note the relation 

where Tg C are the connection coefficients of the n — 1 metric hAB ■ In that space we consider a particle moving 
under the action of the force F % = F u [u, y C )-§^+ F A (u, y c ) g^r- The equations of motion = F l when 
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projected along the direction of u and in the n — 1 space give the equations 

u"-uh AB y' A y' B = F u (74) 

y»A + l u >y>A + T A cy ,B y ,C = F A (?5) 

where u' = 4 s and s is an affine parameter. 

Because the system is autonomous admits the Lie point symmetry <9 t . Using the vector dt and the gradient 
HV H l = ud u we construct two representations of si (2, R) by means of the sets of vectors (see ()36l) .(f4l) )) ) 

d s , 2sd s + ud u , s 2 d t + sud u when fi = (76) 

8 a , -e ±2 ^d s ±e ±2 " s ud u when fx ^ (77) 

and require that the vectors in each set will be Lie symmetries of the system of equations (fT4")l , ([75)1 . In Appendix 
A we show that the requirement of the invariance of the force under both representations ([TBI . (1771) of sl(2, R) 
demands that the force is of the form 

F i = (fu + d u + ±G A (y c ) d A . (78) 

Replacing F l in the system of equations ([73 |) .([75 ]l we find 

u" -uh AB y fA y' B = fu+^G u (79) 
y" A + lu' y ' A +T BC y' B y ' c = ±G A . (80) 

Multiplying the second equation with 2u A h BA y' D and using ((73)) we have 

u 4_^_ (^__„,ID„,IB\ , A „,3 h _ .„.i„,lA„,lD - r>n-„.'D 

ds 

from which follows 



h DB y' D y' B ) + 4u 3 h DA u'y' A y' D = 2G D y lD (81) 

j- (u 4 h DB y' D y' B ) = 2G D y' D . (82) 
ds 

The rhs is a perfect differential if Go = —^.d where is a differentiable function. If this is the case 

we find the first integral 

J = u 4 h DB y' D y' B + 2S (y c ) . (83) 

We note that J involves the arbitrary metric h AB and the function T,(y A ). Furthermore equations I179[) . (|50|) be- 
come 

u" - uh AB y' A y' B = f u +±G u (y c ) (84) 

nA , ^ i , A A /B ic 1 , AB ^ ( C\ / oc \ 

y +~ u y + t bcv y = -—± h ^{y ) tB ■ ( 85 ) 

These are the equations defining the n— dimensional autonomous Riemannian Kepler Ermakov system. 
Using the first integral (|53"|) the equation of motion (1541 is written as follows 

U » = fl i u+ ^M1 (86 ) 

where G = J + G u (y c ) — 2S (y c ) ■ This is the Ermakov-Pinney equation hence we identify ([53")) as the Rie- 
mannian Ermakov integral of the autonomous Riemannian Kepler Ermakov system. 
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5.2 The autonomous Hamiltonian Riemannian Kepler Ermakov system 

In this section we assume that the force is derived from the potential V (u, y ) . Then the equations of motion 
follow from the Lagrangian 

L = \ {u' 2 + u 2 h AB y' A y' B ) - V (u, y c ) . (87) 



The Hamiltonian is 



The equations of motion are 



E = \(u' 2 + u 2 h AB y lA y' B ) +V (u,y c ) . 



u"-uh AB y' A y' M + V, u = o (89) 

y" A + -u'y' A + Y BC y' B y lC + \h AB V, B = 0. (90) 

The demand that Lagrangian (|87[) admits Noether point symmetries which are generated from the gradient HV 
leads to the following cases (see theorem [5] and for details see |12| ) 

Case A: The Lagrangian ([57]) admits the Noether point symmetries ([75]) if the potential is of the form 

V (u,y C ) = ±V (y C ) . (91) 



Ea = ^(u' 2 + u 2 h AB y /A y' B ) + - 2 V( y c ) (92) 





1 


E B 






2 




1 


1+ 






A* 




1 


I_ 






A* 



II 2 1 

.'2 , „,2i, „JA„.IB\ M „,2 | 1 t /■/ / C\ 



The Noether integrals of these Noether point symmetries are 

- 1 (n' 2 +u 2 h AB y' A y' B )+- 2 
h = 2sE-uu' (93) 
I 2 = S 2 E - suu + i« 2 (94) 

where Ea is the Hamiltonian. 

Case B: The Lagrangian ([87)) admits the Noether point symmetries ([77| if the potential is of the form 

V(u,y c ) = -^-u 2 + ^V'(y c ). (95) 
The Noether integrals of these Noether point symmetries are 

„2 _ 1 

(96) 

* s m' + M e 2 "V (97) 
e^ 2 ^ s W + ^e^ 2 ^^ 2 (98) 

where is the Hamiltonian. 

Using the Noether integrals we construct the Riemannian Ermakov invariant Jg, which is common for both 
Case A and Case B, as follows 

Jg = u 4 h DB y' D y' c + 2V (y c ) . (99) 

This coincides with the invariant first integral defined in (|53")l . We note that with the use of the first integral 
the Hamiltonians (jH]) and (JM]) take the form 

E = -u' 2 - ^-u 2 + (100) 
2 2 2u 2 V ' 

which is the Hamiltonian for the Ermakov Pinney equation. 

As it was the case with the Euclidian case of section [5J it can be shown that the Riemannian Ermakov 
invariant is due to a dynamical Noether symmetry [2 2 . 

We collect the results in the following proposition. 
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Proposition 5 In a Riemannian space with metric gij which admits a gradient HV, the equations of motion 
of a Hamiltonian system moving under the action of the potential (fieC) 

V(u,y c ) = -^u 2 + ±V'(y c ) (101) 

admit the sl(2,R) invariance and also an invariant first integral, the Riemannian Ermakov invariant. This 
latter quantity is also possible to be identified as the Noether integral of a dynamical Noether symmetry. 

Without going into details we state the following more general result. 

Proposition 6 Consider a n— dimensional Riemannian space with an r- decomposable metric which in the 
Cartesian coordinates x\, . . . , x r has the general form 

ds 2 = prf-£hdz^dz A + hijdx 1 dx-' , i ) j = r + l,...,n, S = l,...,r (102) 

where r/s\ is a flat non- degenerate metric (of arbitrary signature). If there exists a potential so that the vectors 
e ±fiS ^2a M Dm , are Noether symmetries, where clm are constants, with Noether integrals 

M 

I ± = e ±/» £ a M z' M T fJe^ s ]T a M z M (103) 

M M 

the combined first integral 1 = 1+1- is time independent and it is the result of a dynamical Noether symmetry. 

In the remaining sections we consider applications of the autonomous Riemannian Kepler Ermakov system 
in General Relativity and in Cosmology. 



6 The autonomous Riemannian Kepler Ermakov system in General 
Relativity 

6.1 The autonomous Riemannian Kepler Ermakov system of four degrees of free- 
dom in a FRW spacetime 

Consider the spatially flat FRW spacetime with metric 

ds 2 = du 2 - u 2 (dx 2 + dy 2 + dz 2 ) . (104) 

This metric admits the gradient HV ud u and six non gradient KVs [30] which are the KVs of E 3 . 

We consider the autonomous Riemannian Kepler Ermakov system defined by the Lagrangian (see (jlOll) ) 

(fi.eC) 

L = \ (u 12 - u 2 (x' 2 + y' 2 + z 12 )) + >^u 2 - 1 V (x, y, z) . (105) 
The Euler Lagrange equations are 



" / /•> . >•> . r>\ 2 2V (ar, y, z) 
i + u(x +y + z J — fj u = (106) 



cr// , 2 , V>°(x,y,z) 



x a " + -u'x al - v ' = (107) 



where a = 1,2,3. The Lagrangian (| 105[) has the form of the Lagrangian (|87|l for potential V (u,y 

2 



—u 2 + -^V (x, y, z) hence according to proposition [S] possesses sl(2,R) with Noether symmetries and for 
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both representations (|76[) and (|77|) . The two time independent invariants are the Hamiltonian and the Rieman- 



nian Ermakov invariant (proposition [5]) 

E = l( u ' 2 -u 2 (x' 2 + y ' 2 + z' 2 ))-^-u 2 +^V(x, y ,z) (108) 
2 2 u z 

J Gl = (x' 2 + y' 2 + z' 2 ) + 2V {x, y, z) . (109) 
We remark that if we had considered the representation (|76l) only (that is we had set fi = 0) then we 



would have lost all information concerning the system defined for n ^ 0! We emphasize that in application to 
Physics the major datum is the Lagrangian and not in the equations of motion, therefore one should not make 
mathematical assumptions which restrict the physical generality. 

To assure Liouville integrability we need one more Noether symmetry whose Noether integral is in involution 
with E, Jg 4 . This is possible for certain forms of the potential V (x, y, z) . Using the general results of [TT] where 
all 3d potentials are given which admit extra Noether symmetries we find the following result. 

Proposition 7 The Lagrangian H05\) admits an extra Noether symmetry if and only if the potential V (x, y, z) 
has the form given in table 5, lines: 1,2. table 8, lines :1~4 and table 9, lines:l~4 ofJJljj. 

For example if V (x, y, z) = (x 2 + y 2 + z 2 )™ then the system admits three extra Noether symmetries which 
are the elements of so (3). If V (x, y, z) — Vq then the system admits six extra Noether symmetries (the KVs of 
the three dimensional Euclidian space). 

6.2 The autonomous Riemannian Kepler Ermakov system of three degrees of 
freedom in a three dimensional spacetime 

Consider the three dimensional Lorenzian metric 

ds 2 = du 2 - u 2 (dx 2 + dy 2 ) (110) 

which admits the gradient HV ud u and the three KVs of the Euclidian metric E 2 . In that space consider the 
Lagrangian 

L> = \ (u> 2 - u 2 (x> 2 + y' 2 )) + ^u 2 - ±V (x, y) . (Ill) 

According to proposition [5] this Lagrangian admits as Noether symmetries the elements of si (2, R) . Then from 
proposition [5] we have that the Noether invariants of these symmetries are 

E = l(u> 2 -u 2 (x» + y> 2 ))-£u* + ±V(x,y) (112) 
J G3 = u 4 (x' 2 +y' 2 ) +2V(x,y). (113) 

The requirement that the Lagrangian admits an additional Noether symmetry leads to the condition 
LkvV {x, y) + p = therefore in that case we have a 2d potential and we can use the results of [12] . If 
we demand the new Noether integral to be time independent (p — 0) then the potential V (x, y) and the new 
Noether integrals are given in Table 1 of Appendix B. 

Lagrangians with kinetic term Tk = \ (u /2 — u 2 (x 12 + y' 2 )) appear in cosmological models. In the following 
section we discuss such applications. 

7 The autonomous Riemannian Kepler Ermakov system in cosmol- 
ogy 

In a locally rotational symmetric (LRS) spacetime we consider two cosmological models for dark energy, a scalar 
field cosmology and an f(R) cosmology. 
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7.1 The case of scalar field cosmology 

Consider the Class A LRS spacetime 

ds 2 = -N 2 (t) dt 2 + a 2 (t) e- 2 ^dx + a 2 (t) e ^ (dy 2 + dz 2 ) (114) 

which is assumed to contain a scalar field with exponential potential V (</>) = Vbe~ Cl * , c ^ \/6k and a perfect 
fluid with a stiff equation of state p = p, where p is the pressure and p is the energy density of the fluid. The 
conservation equation for the matter density gives 

p + 6p- = 0^p = ^. (115) 
a a b 

Einstein field equations for the comoving observers u a — jf^dt , u a u a = — 1 follow from the autonomous 



Lagrangian [3"Tl I3"2"] 

We set iV 2 = e c * and the Lagrangian becomes 



L = - 3 "a 2 + ^$ + ^^ 2 - Nka 3 e- C * - N^. (116) 
N 4N H 2N V a 3 V ' 



The Hamiltonian is 



L = (-3d 2 + -a 2 $ 2 + ^-a z 4> z )ae-^ - ka A V e-^ - -f^r. (117) 

4 2 a^e^^ 



E = (-3d 2 + -a 2 {i 2 + k -a 2 ^ 2 )ae-^ + ka 3 V e~^ + -j^- = (118) 

4 2 a d e2 < '' 



If we consider the transformation 

„3 = P x+y ^ _ 

3 V fc 



a 3 = e»-+*, 4> = (x - y) (119) 



where 

x = — !— In ( ^l^-ue*) , y = — !— In ( I±£ue-* J , c = -S= + 1 (120) 

1-c V V2 ) l + c \ y/2 J V6k V ' 

the Lagrangian (|117D becomes 



L = -|m 2 + m ; 

■ ) 



We consider a 2d Ricmannian space with metric defined by the kinematic terms of the Lagrangian, that is: 

ds 2 = ^-6da 2 + ^a 2 df3 2 + ka 2 d<f) 2 ^j ae - ** (122) 

We show easily that this metric admits the gradient HV H l = 6k t c 2 (kad a + cd$) , with gradient func- 
tion H — ^szjfjg fle""^. Therefore the Lagrangian (|117l) defines an autonomous Hamiltonian Riemannian 
Kepler Ermakov system with potential (p ^ 0) 



kVppo 1 
2" ~ ' /i 2 ii 2 



y(u, 2 /' 4 ) = ->V + ^4. (123) 



Because /i ^ this Lagrangian admits sl(2,R) invariance only for the representation (|76p (an additional result 
which shows the necessity for the consideration of the cases p = and p ^ 0!). 
Using proposition [5] we write the Ermakov invariant 

2_, 2 3 , kVapa 1 



J ="VWV + ^' <124) 



1G 



The second invariant is the Hamiltonian E 

P 2. 2 2 (2. a 3 - 2 \ 2 fc^QPo 1 ri9 ,s 

From table 1 line 6 of [IT] we find that the Lagrangian admits three more Noether symmetries 

dp, d z , zd p -/3d z (126) 

with corresponding integrals 

h = u 2 /3 , I 2 = u 2 i , h = u 2 (^jfirzP - ^(3z^j . (127) 
It is easy to show that three of the integrals are in involution, therefore the system is Liouville integrable. 

7.2 The case of / (R) - cosmology 

Consider the modified Einstcin-Hilbert action 



S = J d^x^f—gj (R) (128) 



where f(R) is a smooth function of the curvature scalar R. The resulting field equations for this action in the 
metric variational approach are 33 

f'Rab - \.f9ab + 9a b af - f! ab = (129) 

where /' = ^.y and /" ^ 0. In the LRS spacetime f|l 14[) with N (t) = 1 these equations for comoving observers 
are the Euler-Lagrange equations of the Lagrangian 



L = ( Qaf'a 2 + 6a 2 f"aR - -f'a 3 f3 2 ) + a 3 (f'R - f) . (130) 

The Hamiltonian is: 



E = l^af'a 2 + 6a 2 f"aR - J - a 3 (f'R - f) = 0. (131) 

Again we consider the 3d Riemannian space whose metric is defined by the kinematic part of the Lagrangian 

(pro 

ds 2 = Uaf'da 2 + 12a 2 f" da dR - 3a 3 f'd(3 2 . (132) 
This metric admits the gradient HV 

H" = \ Ud a + j^dj^J (133) 

with gradient function H — 3a 3 /'. 

In order to determine the function f(R) we demand the geometric condition that Lagrangian (|130|) admits 
s(2, R) invariance via Noether symmetries (see [H]). Then for each representation ([TB")) . (|77jl we have a different 
function f(R) hence a different physical theory. 

The representation (|76[) in the present context is: 

1 ( , f' a \ ,2, ' ( ., , /' 



d t , 2td t + - \ad a + jjjd R j , t 2 d t + - (^5 Q + jjjd R j . (134) 
The Noether conditions become 

- 4a 3 f'R + 7 -a 3 f + p = 0. (135) 
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These vectors are Noether symmetries if p — and 

f(R) = R Z s. (136) 

however power law / (R) theories are not cosmologically viable [33] . 

The second representation (1771) in the present context gives the vectors 



8 t , -^d t ± \e ±2 ^ (ad a + j^d R j . (137) 
The Noether conditions give 

-4a 3 f'R+ ^a 3 f + 3fi 2 a 3 f +p = 0. (138) 

These vectors are Noether symmetries if the constant p = and the function 

f(R) = (R-2A)i (139) 

where 2 A = 3/i 2 . This model is the viable A;, c CDM-like cosmological with b = l,c = |. |35j . 

We note that if we had not considered the latter representation then we would loose this interesting [§ result. 
For the function (|139l) the Lagrangian (|130[) becomes for both cases (if A = we have the power-law 

f(R)=Ri) 

21 i o 21 o s • 

L = — a (R- 2A)~* a 2 - —a 2 (R- 2A)~* aR 

8 8 (i?-2A) 5 



Furthermore there exist a coordinate transformation for which the metric (|132l) is written in the form of (187 
We introduce new variables u, v, w with the relations 



In the new variables the Lagrangian (1140p takes the form 



21 \ 3 e 12v 

ueF , R = 2A+ — -, $ = V2w. (141) 



L=iu 2 -iu 2 (u 2 + ™ 2 ) + ^u 2 ~- — . (142) 
2 2 V ; 2 42 u 2 V ' 



The Hamiltonian (|13ip in the new coordinates is 

2 2 V ; 2 42 u 2 V ' 

The Lagrangian (1142)) defines a Hamiltonian Riemannian Kepler Ermakov system with potential 

from which follows the potential V(v) = j^e 12 ". In addition to the Hamiltonian the dynamical system admits 
the Riemannian Ermakov invariant 

j f=u ^{f + w 2 ) + ^e 12 \ (144) 

The Lagrangian (I142p admits the extra Noether point symmetry d w with Noether integral I w = u 2 w (see 
Table 1). The three integrals E, I w and J/ are in involution and independent, therefore the system is integrable. 



3 The importance of the result is due to the fact that it follows from a geometric assumption which is beyond and above the 
physical considerations. Furthermore the assumption of Noether symmetries provides the Noether integrals which allow for an 
analytic solution of the model. 



18 



8 Conclusion 



In this work we have considered the generalization of the autonomous Kepler Ermakov dynamical system in the 
spirit of Leach [5], that is using invariance wrt the sl(2,R) Lie and Noether algebra. We have generalized the 
autonomous Newtonian Hamiltonian Kepler Ermakov system to three dimensions using Noether rather than 
Lie symmetries and have determined all such systems which are Liouville integrable via Noether symmetries. 
We introduced the autonomous Riemannian Kepler Ermakov system in a Riemannian space which admits a 
gradient HV. This system is the generalization of the autonomous Euclidian Kepler Ermakov system and opens 
new fields of applications for the autonomous Kepler Ermakov system, especially in relativistic Physics. Indeed 
we have determined the autonomous Riemannian Kepler Ermakov system in a spatially flat FRW spacetime 
which admits a gradient HV. As a further application we have considered two types of cosmological models 
which are described by the autonomous Riemannian Kepler Ermakov system, i.e. scalar field cosmology with 
exponential potential and f(R)— gravity in an LRS spacetime. 
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Appendix A: Proof of the force in equations (1741). (1751) 



We require that the force admits two Lie symmetries which are due to the gradient HV H = ud u (if we require 
the force to be invariant under three Lie symmetries which are due to the gradient HV then it is reduced to the 
isotropic oscillator). From theorem Q] (for details see [12]) we have the following cases. 

(I) Case fi = 

In this case the Lie symmetries are 

d s , 2sd s + ud u , s 2 d s + sud u . 
The condition which the force must satisfy is 

L H F l + dF l = 0. 

Replacing components we find the equations 

(*-!)*«■ = 

^-F A )u + dF A = 
ou J 

from which follows 

p u — \ p u p A — —F A 

u (d-l) ' u d 

Because the HV is gradient, condition A2 of theorem [1] applies and gives the condition 

L H F l +AF { + a 1 H i = 
from which follows a\ — and d = 4. Therefore 

F u = l_ G u {C) F A = ± G A {C) 

(II) Case pL ^ 

In this case the Lie symmetries are 

-e ±2 ^d s ±e ±2 ^ud u 
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The condition which the force must satisfy is 

L H F i + AF l + a\H % = 
We demand a± ^ and obtain the system of equations: 

| -0-F u ) u + 3F u + am = 
\ou J 

( u + 4F a = 
\ou J 

whose solution is 

F u = f i 2 u + \g u , F A = — A G A 

where we have set a\ = — 4/i 2 . 

Appendix B: Table 1 



Table 1: Potentials for which the Lagrangian (jllip admits extra Noether point symmetry 



Noether Symmetry 


V(x,y) 


Noether Integral 


d x 


f(y) 


I x = u 2 x' 


8y 


f{x) 


I y = u 2 y' 


yd x - xdy 


f(x 2 + y 2 ) 


I xy = u 2 (yx 1 - xy') 


d x + bd y 


f{y- bx ) 


hby = u 2 (x' + by') 


(a + y ) d x + (b - x) d y 


f (i (x 2 + y 2 ) + ay - bx) 


Iabxy = {a + y) u 2 x' + (b-x) u 2 y' 


d x , d y , yd x - xdy 


V 


Ix i Iy i -^xy 
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